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Abstract 

We perform a multi-parameter global analysis of all data available till date from the 
ATLAS, CMS and Tevatron experiments, on the signals of a Higgs boson, to investigate 
how much scope exists for departure from the standard model prediction. We adopt a 
very general and model-independent scenario, where separate deviations from standard 
model values are possible for couplings of the observed scalar with up-and down-type 
fermions, W-and Z-boson pairs, as well as gluon and photon pair effective interactions. 
An arbitrary phase in the coupling with the top-pair, and the provision for an invisible 
decay width for the scalar are also introduced. After performing a global fit with seven 
parameters, we find that their values at 95% confidence level can be considerably 
different from standard model expectations. Moreover, rather striking implications of 
the phase in top-quark coupling are noticed. We also note that the invisible branching 
ratio can be sizeable, especially when the couplings of the Higgs to W-and Z-pairs are 
allowed to be different. 
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1 Introduction 



The recent announcements on Higgs boson search from the 2012 run of the Large Hadron 
CoUider (LHC) [H E], in conjunction with earher resuhs, have brought a refreshing gale 
into the stuhified arena of particle physics, starved of breakthrough for several decades. 
Physicists tend to concur that both the CMS and ATLAS collaborations have found a spin- 
zero bosonic resonance of mass 125 - 126 GeV. However, the issue of whether it is exactly the 
Higgs boson predicted by the standard electroweak model (SM) [3lll] is still somewhat open. 
This question can be settled by the close examination of accumulating data over a longer 
period, thus revealing not only the complete picture of the electroweak symmetry breaking 
sector but also any secret message of new physics buried within the available results. 

Any departure from SM predictions in the observed scalar is ultimately reflected in its 
interactions with pairs of fermions or weak gauge bosons, and also the loop-induced effective 
couplings to photon and gluon pairs. However, these interactions are rather intricately 
twined in the calculation of rates for various final states into which the Higgs can decay. The 
trace of non-standard physics there can be looked for in two ways — either in the context of 
specific theories where the different new couplings are all dictated by the model parameters, 
or in a phenomenological, model-independent analysis of as general a nature as possible. 
We attempt to add the present investigation in the second direction, to the already growing 
volume of extant studies in the light of the data piling up |5l El [7] . 

We start by assuming that the SU{2)l x U{1)y gauge invariance holds at the energy 
scale under scrutiny. Under such circumstances, the aforementioned interactions of the 
Higgs boson can be different from the SM expectations due to (a) the observed scalar having 
admixtures with other states, or (b) the presence of additional, gauge invariant effective 
operators which contribute to the same couplings. In either case, a renormalisation of the 
SM couplings may take place, and interactions with new Lorentz structure may also appear. 
Only the former possibility is addressed in the present work. 

The modification of the SM couplings of the Higgs is subject to various experimental 
constraints, the severest of them often coming from electroweak precision observables. Our 
purpose is to investigate, within such constraints, how much allowance for departure of 
various interaction strengths from the SM values can be made by the currently observed 
rates in different final states, as found from a global fit of these strengths based on chi- 
square minimisation. The available best fits for the rates in different channels, scaled by 
the corresponding SM expectations, are our inputs, and we use them to obtain the widest 
allowed ranges for various modified couplings as well as an invisible decay width, with the 
hope that this will guide us to the direction where new physics may lie. 

We start by assuming that the Higgs (or, to be more precise, the observed boson, to 
which we refer here as the Higgs) can differ from the SM expectations in all respects — 
couplings with T3 = -1-1/2 and T3 = —1/2 fermions (the departures being different in the 
two cases), couplings with both W-and Z-boson pairs (again with potentially different devi- 
ations), and also the effective couplings to photon and gluon pairs, where additional effects 
over and above the modified fermion and boson couplings are included as possibilities. While 
the modifications mentioned above can be described in terms of a given framework (such 
as a chiral Lagrangian), we deliberately take them as completely free and independent pa- 
rameters, to make our study free of any theoretical bias. Besides, it is also assumed at the 



2 



beginning that the Higgs can have a finite invisible decay width. We use the CMS as well 
as ATLAS results in various channels for both the 2011 (7 TeV) and 2012 (8 TeV) runs [S]. 
Results from the Tevatron, wherever available, are also used as data points in obtaining the 
least-square fits |9j. 

As has been already mentioned, a number of similar investigations have appeared in the 
literature [7]. Though they are all instructive, the present study may be of particular use in 
the following respects: 

• We start by taking all of the couplings, tree-level as well as loop-induced, to be unre- 
lated and free parameters. The couplings to up-and down-type quarks, as also those 
to W-and Z-pairs, are allowed to be different at the same time. Most earlier studies 
have not allowed uncorrelated variation of HWW and HZZ interactions, albeit in a 
gauge-invariant fashion; also, some of them set the interactions with up-and down-type 
quarks to identical value. 

• We also investigate the effect of an arbitrary phase in the fermion pair couplings, 
relative to the W-pair couplings, which is also varied as a free parameter in obtaining 
the global fits. 

• We do not restrict ourselves by assuming the anomalous interactions to arise from some 
underlying scenario like a chiral Lagrangian, which implies concomitant variation of 
different couplings. 

• We take into account the possibility that the effective Hgg and H'yy couplings are 
modified due to effects other than non-standard interactions of the Higgs with W- 
and fermion pairs. Furthermore, the possibility of loop-induced Higgs decays being 
modified differently due to coloured and colourless fermions running in the loops is 
taken into consideration. We keep the modifications due to these two different effects 
separate and uncorrelated. 

• An invisible decay width for the Higgs is allowed as a free parameter. 

• We parametrize all new physics effects at the coupling level (excepting for the in- 
visible width, where we take the width as a free parameter, for the sake of model- 
independence). This is in contrast to cases where branching ratios of the Higgs into 
some channels are taken as free parameters, which results in various new physics ef- 
fects getting entangled, since branching ratios also involve the total decay width of the 
Higgs, where all (modified) interactions contribute. 

• In obtaining the 2a ranges of allowed values for various couplings, we have fixed the 
remaining parameters not at their SM values but at their global best fits corresponding 
to the minimum of the function. This makes our analysis completely unbiased, from 
the angle of physics beyond the SM. 

• Even in channels where the gluon fusion channel for Higgs production dominates, the 
vector boson fusion (VBF) and associated production channels contaminate the rates 
to varying degrees. Also, the VBF and associated channels have both HWW-a.nd 
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HZZ-induced contributions. These warrant a careful treatment if one is allowing 
uncorrelated deviations from SM in the HWW and HZZ interactions. Difficulties on 
this front are compounded by different efficiency of cuts in different channels. Such 
issues have been taken into account in our analysis. 

• In cases where only the 2011 data are used, where the best fit values for various 
rates come with asymmetric error-bars, the asymmetry is retained in the analysis that 
follows. 

In section 2, we outline our parametrization of the various new physics effects mentioned 
above, and discuss the motivations of their origin. The details of input data sets from the 
LHC and the Tevatron used in our global numerical analysis, and the methodology adopted 
have been described in section 3. We discuss our results of the best-fit values obtained and 
the allowed confidence intervals for the various parameters in section 4. We summarise and 
conclude in section 5. 



2 New physics effects: parametrization 

In this section, we describe the parameters used to encapsulate any likely new physics effect 
hidden within the data on the Higgs. We also outhne the motivations for choosing these 
parameters, and indicate the ranges over which we let them vary to obtain their favoured 
locations in the light of observations till date. 

Fermion couplings 

Classifying all T3 = -1-1/2 fermions as u, and all T3 = — 1/2 fermions as d, we assume 

r iS TT- 

V 

^Hdd = aa—Hdd (1) 

V 

where H denotes the scalar with a Higgs-like appearance. This parametrization implicitly 
assumes the couplings of H to all fermions to be proportional to their masses. The assumed 

difference in interaction strengths to up-and down-type fermions include the possibility of H 
being part of a scenario containing more than one doublets, including the supersymmetric 
case. In general, a relative phase between the couplings is assumed. This phase allows 
various degrees of interference among fermion loops, and between the fermion and W-boson 
loops, for example, in the decay H — )■ 77. As will be seen in the next section, we perform 
a scan over the phase 5 like we do over au and ad- Even in cases where we neglect a non- 
trivial phase, we allow a;„ and Ud to be both positive and negative, in order to account for 
constructive as well as destructive interferences. Moreover, the phase 5 enters seriously into 
Higgs phenomenology only via the top quark couplings. Therefore, nothing is affected by 
dropping it for the first two families, in case it is subject to any constraints from flavour 
physics. 
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There are essentially no limits on and \ad\ from earlier data, except those from 
perturbativity of the Yukawa couplings. Keeping this in mind, we take the maximum value 
of la^l to be 2 (keeping the top quark Yukawa coupling in mind), while la^l is prima facie 
allowed to lie all the way upto 40. However, even before the full analysis is done, we find 
that the rate for H — > 77 for large \ad\ gets suppressed well below limits admissible by 
current observations, as the bb mode then dominates overwhelmingly. With this in view, 
we also limit the maximum value of {a^l to 2 in our final analysis. Since this renders the 
contribution of all T3 = —1/2 SM fermions to r{H — 77) negligibly small, we also do not 
miss anything by dropping the phase in ad- 

Gauge boson pair couplings 

We parametrize the interactions of the observed scalar to a pair of weak gauge bosons as 



V ^ 
2 

Chzz = Pz—HZ,Z^ (2) 

V 

where the Lorentz structure of the interaction has been tentatively taken to be the same as in 
the standard model. It should be emphasized that this parametrization, especially allowing 
Pw 7^ Pz, allows one to address the apparent suppression of the WW channel as opposed 
to the ZZ channel, as evinced in the 2011 data. Such anomalous interactions can arise if, 
again, the Higgs has admixtures of other doublets or scalars in some other representations 
of SU(2), and also via loop effects in specific models. 

Clearly, one faces precision electroweak constraints here, in particular, for [3w 7^ [3z, 
which entails a breakdown of custodial SU(2), and is thus restricted by the T-parameter. 
Such anomalous couplings can arise, for example, from gauge invariant effective operators, 
an example being [10] 

0<^ = ^(D,<|.)t$<|.t(Z}^$) (3) 

where <l> is a set of SU(2) doublet scalars, out of which H is the lightest mass eigenstate, 
and A is the scale below which the effective operator is defined. This operator in itself gives 
rise to unequal Pw and Pz- However, taking this operator alone, precision constraints yield 
the hmits [TOl [II] 



0.991 <(5w^ 1.001 (4) 
0.997 ^(5z^ 1.028 (5) 

It is hardly expected to see any appreciable effects of such variation on observable rates. 
However, one can have less constrained couplings if one includes other effective operators 
which, however, give rise to additional HWW and HZZ interactions involving derivatives. 
We do not rule out such possibilities, but neglect the effects of derivative couplings for the 
time being, and vary (3w and Pz between and 2, in a purely phenomenological way. An 
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analysis including the derivative couplings will be presented in a subsequent study. We also 
consider the case where /3w = Pz = thereby restoring tree-level custodial invariance. In 
this case, for a Higgs mass of 125 GeV, electroweak precision constraints restrict /3 in the 
range |5j 

0.84 < < 1.4 (6) 



Effective gluon-gluon and photon-photon couplings 

The gluon fusion channel is the dominant production mode for a Higgs of mass around 
125 GeV, and is overwhelmingly driven by the top quark loop. Therefore, a departure of 
Q!„ from unity will be reflected in the production cross-section (though the phase 5 will be 
ineffectual). Similarly, the two-photon amplitude, which is responsible for perhaps the most 
discussed final state nowadays, has contributions from both fermion-and VT-induced loops. 
Thus the parameter also dictates the rate for the two-photon final state. 

This is, however, not the entire story. Both of the aforementioned loop-induced pro- 
cesses can have modified contributions, beyond the coverage of the «-and /3-parameters, if 
additional states contribute in the loops. The most obvious example is the contribution 
of Kaluza-Klein towers in theories with extra compact dimensions, where fermions and/or 
gauge bosons propagate in the bulk. Due to such (and perhaps other) possibilities, it is 
necessary in a general study to include an additional parameter to properly quantify new 
physics effects in the gluon fusion channel. For the two-photon amplitude, this parameter can 
well be different, since new physics can be quite different in the coloured and non-coloured 
sectors. 

With this in view, we parametrize the gluon-gluon-Higgs and Higgs-photon-photon am- 
plitudes as follows: 



L^^ = -x^g{au, ad, /3w, 5)-^HFf,^F^"' (7) 

where Xg and encapsulate the overall modification due to new intermediate states in the 
two cases. The functions /(a^) and g^au^cxd, Pw,^) encapsulate the modifications of these 
couplings due to fermion and W-boson loops. We shall discuss their detailed forms in the 
next section when we discuss the departures of the Higgs production and decay widths from 
their SM values in detail. Since there is no restriction till now on Xg and x^ , we let each of 
them vary from 0.2 to 3.0. The lower and upper bounds have been set keeping in mind that 
we do not have too much room for varying them beyond a range, and still being consistent 
with the data. 

Invisible width 

Since the earlier analyses of the 2011 data have led to different conclusions about a 
possible invisible width of the Higgs |6j, we keep this possibility alive in our analysis. Such 
an invisible width can occur if, for example, the Higgs serves as a portal to a 'dark matter' 
sector. The exact expression of the width in terms of the coupling to the dark sector will 
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require one to know the nature of the invisible particle(s), for example, whether they are 
scalars or spin-1/2 objects. In order to be model-independent, we take as a free parameter 
the invisible width Tmv, which is independent of the nature of the invisible state, and is also 
not entangled with other new physics effects. 

Since there is very little guideline on the range over which Ymv should be varied in order 
to obtain the value corresponding to the minimum of chi-squared, we start from e, the 
invisible branching ratio, and let it vary between and 1. In each case, the invisible width 
is expressed as 



where Yl ^vis is the total decay width into all visible channels. 



3 Methodology of analysis 
3.1 Input data 

Table 1 contains the details of all the data points used in our analysis. This includes the 
combination of 7 TeV (with 5.1 fb~^ data) and 8 TeV (with 5.3 fb~^ data) results from 
CMS [HIH] in all channels, namely, 77 (inclusive), ZZ* Ai,WW* ££z/z/, 777 j, r+r~ and 
bb. For the ATLAS experiment |2l [8], the 77 and ZZ* results are available as similar 7 + 8 
TeV combinations, whereas the 2011 results alone have been used for the remaining channels. 
The integrated luminosities for ATLAS are 4.9 fb"^ and 5.9 fb"^ for the 7 and 8 TeV runs 
respectively. The 'j'jjj results from ATLAS are not yet available. Furthermore, we have 
used the Tevatron (combined CDF and DO) results for WW*,'~f'~f and bb, for an integrated 
luminosity of 10fb~^ f9j. Thus we have fourteen input data points altogether for our global 
analysis. All the SM production cross-sections and decay widths for the Higgs have been 
taken from the results reported by the LHC Higgs Cross Section Working Group |12j . 

In calculating the modifications of various branching ratios of the Higgs, we have used 
mt = 173.5 GeV, = 4.7 GeV, = 1.777 GeV and mw = 80.4 GeV \13\. The best fit 
value for the Higgs mass {rriH) reported by CMS is 125.3 ± 0.6 GeV P, whereas it is 126.5 
GeV [2] for ATLAS. The average of these two values, namely, 125.9 GeV has been used by 
us, in line with, for example, Ref. [H]. For the Tevatron analysis, tuh = 125 GeV has been 
used, following combined CDF and DO analysis as reported in Ref. [9]. 



Experimental collaborations have reported various observed signal strengths in the i chan- 
nel in terms of fii = a"^'^/af^^, with cxj as the respective uncertainty. Here, cr°^* refers to 
the observed signal cross-section for a particular Higgs mass, while af is the signal cross- 
section for an SM Higgs with the same mass. We calculate the corresponding values of yUj 
for various points in the space spanned by the parameters in terms of which we have tried 




(8) 



3.2 Methodology 
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Experiment 


Channel 




Mh (GeV) 


Energy (TeV) 




H W+W- 


"-'•'-'^-0.32 






Tevatron 




1 07+0.74 
-0.68 


125 


1.96 




if — >■ 77 


3.6211^^ 








if — >• 77 


1 q+0.50 

-■-•^-0.50 


126.5 


7 + 8 




H ZZ* Al 


1 q+0.60 
-'-•'-'-0.60 


126.5 


7 + 8 


ATLAS 


H WW* l+l-uP 


"-'•'-'^-0.60 


126 


7 




H ^ Tf 




126 


7 




H^bb 




126 


7 




ii" —)■ 77 


1.00_Q 43 








H ZZ* Al 


n 7+0.40 
' -0.40 








H WW* l+l-vV 








CMS 


H ^ Tf 


-O.lAtrrt 


125.3 


7 + 8 




H^bb 


'-'■^^-0.66 








H T133 


-'-•'J°-1.06 







Table 1: Input data set used in our analysis, with the values of /2j in various channels and 
their la uncertainties as reported by the ATLAS, CMS and Tevatron collaborations. 

to capture the departure of the Higgs interactions from their SM values, as explained in the 
previous section. One can express /Xj as 

l^i = i^f"^ X Rf<^v/E"idth (9) 

where RY"'^-, R'^'^"''^ and Ry'^'^^^ are the factors modifying the corresponding SM production 
cross-sections, decay width in a particular channel, and the total decay width of the Higgs 
respectively. 

The relevant production mechanisms at the LHC and Tevatron for the various channels 
are listed below: 

• For 77 (inclusive), ZZ* — )■ Ai,WW* — )■ iiuu and r+r^, one has to include all the 
production processes, namely, gluon fusion, vector boson fusion, associated WH and 
ZH production, and Higgsstrahlung from top-antitop pairs. 

• Only the associated WH and ZH production channels can lead to bb final states that 
can be separated from backgrounds. 

• VBF and GF are the dominant production modes for the 77jj channel reported by 
CMS; in particular, the contribution of VBF dominates once the appropriate cuts 
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to identify forward tagged jets are imposed. This is reflected in the corresponding 
efficiencies: 15% for VBF, and 0.5% for GF. 



In the most general case, the production cross-sections in the gluon fusion (GF), as- 
sociated production with a Z (ZH), associated production with (WH) and associated 
production with ti (tiH) are modified by the factors x^a^ (Rgf), /3| (Rzh), Pw {Rzh) and 
«u [Rttn) respectively. In the vector boson fusion (VBF) channel, the corresponding factor 
is given by 

RVBF ^ ^Jk±^ (10) 

The factor Rvbf requires some explanation. In order to obtain this, first of all, we have 
used the fact that the interference of the WW-fusion and the ZZ-fusion diagrams is of the 
order of 1%, and can therefore be ignored in our calculation. Secondly, the WW-fusion 
contribution to the total cross-section is roughly 3 times that of the ZZ-fusion contribution. 
For details on this point we refer the reader to Ref. [15]. We have cross-checked these facts 
for the LHC energies using the VBF@NNLO code of Bolzoni et al [16]. 

Note that these factors can only be used for inclusive channels of Higgs search, whereby 
the efficiencies for event selection are the same in various production modes. For channels in 
which special kinematic selection criteria are used, one has to also include the corresponding 
efficiency factors, as explained below for the 77jj channel included by CMS. 

For the 77 jj channel reported by CMS, we have considered VBF and GF as the two 
most dominant production modes [TT] . As reported by CMS, the overall acceptance times 
selection efficiency of the dijet tag for Higgs boson events is 15% for VBF {^vbf) and 5% 
for GF {^gf)- Therefore, for the VBF channel, the factor modifying the SM production 
cross- sect ion is given by 

p _ ^VBfRvBFCvbF + ^GfRgF<^gF (-,-,\ 
WbfO^VBF + ^GF<^GF 

Next, we consider the Higgs decay widths in the channels under study. The Higgs decay 
widths in the ZZ* , WW*, rf, bb, cc and gg channels get multiplied by /3|, f3^, a^, and 
x^gO^ respectively. In the loop-induced 77 channel, the contribution due to the top quark 
and W-boson loops (and a small contribution due to the bottom and tau loops) to the Higgs 
decay width is modified by the factor given by 

^ ||a„e-^Af^,(n) + ladA'^/^in) + adA'^/^jr,) + PwA^^w)? 
"'"^^ l|^f/2(n) + Mf/,(r.) + Af/,(r.) + Af(r^)P ^ ^ 

where the relevant loop functions are given by [15] 



A^,,{n) = 2[n + {n-l)f{n)]rr^ 

A^{n) = -[2rf + 3r, + 3(2r, - l)/(r,)]r-2 (13) 



Here, /(tj), for Tj < 1 is expressed as. 
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while, for > 1, it is given by 



(14) 




1+ y/T^ 

1- 



-1 



1 



— m 



(15) 



In the above equations Tj denotes the ratio m^j/Amf. 

With the values of thus calculated for various points in the parameter space, we first 
obtain the best fit values for these parameters (upto 7 at a time) corresponding to the global 
minimum of the function x^, defined as 



Note that, the experimental collaborations, in some cases, have reported asymmetric 
error bars on the data. In order to include such error bars in the above definition, we use the 
following prescription. If (/Xj — /2j) > 0, we use the positive error bar af, while if {ni — fii) < 0, 
we use the negative error bar a~ [IH]. Note that, in cases for which we have combined more 
than one experimental data points to obtain a single input data, we obtained the average 
signal strength fi and the corresponding uncertainty a using the following relations: 



This method has been used, for example, in case of CMS data on inclusive 77 and 77jj 
channels, where the results were presented in four different categories in the former case 
while with both loose and tight dijet tags in the latter case [T7] . 

We have also used the above method in combining different contributions to the theoret- 
ically calculated ^ values, for example, while combining the contributions of WH — )■ Ivhb, 
ZH — )■ l~^l'~bb and ZH — )■ uubb to associated Higgs production with gauge bosons and the 
subsequent decay of the Higgs to a bottom pair. Such combinations were necessitated by 
the fact that the experimental collaborations have reported a single signal strength value in 
the bb channel. 

After obtaining the best-fit values for the parameters by minimizing the function, we 
consider two-dimensional contours for various pairs of them about the global minimum, fixing 
the remaining ones at their best-fit values. The contours are drawn for 68.3% and 95.4% 
confidence intervals. One important point in which our study differs from most earlier ones is 
that, in the cases where these two-dimensional contours are drawn, the remaining parameters 
are fixed not at their SM values but at values corresponding to the global minimum of the 
fit. The various values of Ax^ for the 68.3% and 95.4% confidence intervals are tabulated, 
for example, in Refs. [191 120] . 




(16) 





(17) 
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Case 


Pw 


Pz 






Xg 




e 


6 


A 


1.2 


1.6 


-0.6 


-1.2 


1.8 


1.3 


0.6 


0* 


B 


1.06 


1.06 


-1.05 


0.95 


0.8 


1.0 


0.2 


0.55 


C 


1.07 


1.07 


-1.3 


-0.96 


0.65 


0.95 


0.2 


0* 



Table 2: Best-fit values of the various parameters in the three cases considered. In cases A 
and C, 6 has been fixed at (indicated with a '*'). In cases B and C, the relation = Pz 
has been imposed, and their values have been restricted within precision constraints. 

4 Results 

We obtain the global fits in three different cases. In the first case (case A), the phase 6 is 
set to zero, and the seven remaining parameters are varied across the pre-decided ranges, as 
described in section 2. In the second case (case B), we try to see the effects of the phase, 
and set (3w = /3z in the process. The same exercise is done in the third case, but with the 
phase set to zero again. The simplification in the last two cases is done because f3w Pz 
amounts to the breaking of custodial SU(2) at the tree-level, and is subject to rather stringent 
constraints from precision electroweak observables. Though we go beyond these constraints 
in the first case, assuming the simultaneous existence of more than one gauge-invariant 
effective interactions at the same time, we desist from such speculative exercise in the cases 
B and C, and keep Pw = Pz within the most stringent precision electroweak limits. It should 
be noted that the phase 5 does not affect the loop-induced contribution to the T-parameter, 
since it cancels in the relevant self-energy diagrams. Table 2 contains the best fit values of 
various parameters in the three cases. 

The obtained best-fit values for the parameters have been tabulated in Table [2J The 
minimum values of x^, Xmm 6.19, 7.03 and 7.03 for cases A, B and C respectively. Note 
that, the number of degrees of freedom in our global fits is 7, since we have 14 input data 
points, and 7 parameters. We present the two-dimensional 68.3% and 95.4% confidence 
interval contours for various pairs of parameters about the global minimum, for cases A and 
B, in Figures [1] and [2] respectively. As mentioned before, while drawing these contours, we 
have fixed the remaining parameters at their best-fit values. In this connection, we note that 
the consistency of our methodology has been checked by also performing a two parameter 
fit of the data (with the parameters pertaining to the overall modifications of the fermion 
and gauge couplings), and our results in that case agree with those presented in recent 
hterature [7]. 

Let us first consider the case with 5 = 0. First of all, there is a global tendency in 
the experimental results reported so far towards best fit values of the WW rates that are 
consistent with a shortfall compared to the SM prediction. For ZZ, on the other hand, there is 
an excess for ATLAS and a shortfall, but less than that for WW, for CMS. Though these are 
all consistent at the la level, the least-square minimisation inevitably yields a best fit for Pz 
higher than that of Pw Since this causes one to go beyond the most stringent of precision 
constraints, one has to admit that such a best fit central value favours the simultaneous 
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Figure 1: Two-dimensional contour plots for 68% and 95% confidence intervals, for case A, 
with rest of the parameters fixed at their best-fit values. The best-fit point is also marked 
separately by a In this case 6 has been fixed at 0, whereas < /3w,/3z < 2.0, and 



presence of more than one gauge-invariant anomalous WW/ZZ interactions for the Higgs. It 



12 




Figure 2: Two-dimensional contour plots for 68% and 95% confidence intervals, for case 
B, with rest of the parameters fixed at their best-fit values. The best-fit point is also 
marked separately by a In this case 6 has been varied in the range {0,7r}, whereas 
0.92 <I3< 1.18, with (3 = (3w = Pz- 



should, of course, be kept in mind that one has to wait for the accumulation of more data 
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before a final verdict can be spelt on this. In cases B and C, where we set 0.92 < /3 < 1.18, 
with (3 = Pw = /Sz, we find the best-fit values of (3 to be close to 1. 

As far as the best fit values for the parameters multiplying the fermion couplings are 
concerned, we find that receives a sign relative to [3wi with a value —0.6, —1.05 and —1.3 
for cases A, B and C respectively. This relative sign between the top-quark Yukawa coupling 
and the W-boson coupling to the Higgs can be traced to the fact that the observed excess in 
the 77 mode by both the ATLAS and CMS experiments can be explained by a constructive 
interference between the loop amplitudes due to the top and the W-loop. In contrast to 
the magnitude of the best-fit value of ad in all cases turns out to be close to unity. Its sign 
does not affect the Higgs production cross-sections and decay branching ratios much, except 
for a small effect in the 77 decay mode. However, to be completely general, we varied it 
in the range —2.0 to 2.0. One should note that, as can be seen in Figure [H the 68% and 
95% contours allow for both positive and negative values of a^, and these contours bring 
out the global picture more clearly than just the best-fit values. It should be noted that ad 
has a positive best fit value only when the phase 5 is turned on. This elicits an interesting 
feature of the global fit, when the phase in fermion couplings is included. We also note that, 
from contours involving au and ad in Figure [T] and in Figure [H it is clear that values 
close to zero of these two parameters is clearly disfavoured. This, therefore, should severely 
constrain fermiophobic models for the Higgs j21] . 

We also find the best-fit values for Xg and to differ, and at the same time lie away 
from their SM value of 1.0. This suggests that not only is the presence of new coloured or 
electrically charged states coupling to the Higgs plausible, but that the number of states 
contributing to the gluon fusion process can be different from those contributing to the 77 
decay mode. 




Figure 3: Variation of the function with the invisible branching fraction of H (e) in cases 
A (left panel) and B (right panel). In case A, 5 = and /3vk 7^ Pz-, whereas in case B, 5 has 
been varied in the range {0, vr} and 0.92 < /3 < 1.18, with {3 = f3w = (3z- 



In case A, we find a rather large allowed value for the invisible branching ratio of the Higgs 
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boson. This may be due to the fact that the production rate in the gluon fusion channel, 
as well as the Z boson contribution to the associated production and VBF processes (as the 
best-fit of fiz is 1.6) are enhanced. This creates room for the simultaneous enhancement of 
certain visible modes (77, ZZ*^ as well as the existence of substantial invisible branching 
fraction. Note that we observe a reduction in e once fiw and fiz are set to be equal, and in 
the range allowed by electroweak precision data (cases B and C), the best-fit value in both 
these cases being 0.2. Thus our conclusion is that, while invisible decays of the Higgs are 
more constrained once a correlation between the HWW and HZZ couplings is enforced, the 
overall freedom in a seven-parameter space nonetheless keeps enough scope for such decays. 
We demonstrate the variation of the function with e, in cases A and B, with rest of the 
parameters fixed at their best-fit values, in Figure [31 




Figure 4: Variation of the function with the phase in the up-type quark Yukawa coupling, 
6, in case B. In this case 6 has been varied in the range {0,7r}, whereas 0.92 < f3 < 1.18, 
with (3 = (3w = (3z- 



The interesting feature to be noted is that the phase 6 is allowed to assume values as large 
as 0.55 in the best fit. Once this non-trivial phase is there, all other anomalous interactions 
tend to come close to the SM expectations, the only exceptions being and e. However, we 
draw the attention of the reader in this context to Figure |H Although the best fit value is 
0.55, the function is very slowly varying in the entire range between and 0.6. The curve 
as a whole indicates that the best fit value is relatively unspecified between and 1, and we 
have to wait for more data to have a clearer picture about the occurrence of this phase. 



5 Summary and conclusions 

The attempt in this study is to see how much scope of new physics is contained within the 
data on Higgs search, taking into account not only the results available from 7 and 8 TeV 
runs of the LHC but also the full data set of the Tevatron. We have taken a completely 
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model-independent stand, without any bias of correlation between the anomalous Higgs 
couplings of T3 = +1/2 and —1/2 fermions, as well as the W-and Z-couplings. We have 
in addition accounted for additional states contributing to the effective interaction of the 
Higgs with photon and gluon pairs. Since such contributions can come through different sets 
of states for the gg and 77 pairs, we assign a separate uncorrelated parameter for each of 
these processes. Furthermore, we have taken into account an arbitrary phase in the top-pair 
coupling with respect to that from a W-pair, which can in principle non-trivially affect the 
loop-induced decay H — > 77. Last but not the least, a non-vanishing width for the Higgs 
decaying into invisible final states is kept as a free parameter in our global fit. 

Our study takes into account not only the contributions to all final states from the dom- 
inant gluon fusion channel but also the VBF and all associated production subprocesses, 
with appropriate weightage to the efficiencies for different channels. Also, asymmetric un- 
certainties in the data, as reported by the experiments in some cases, have been used in our 
analysis. 

The global fits for the minimum of over seven-parameter scans, for three sets of 
combinations, yield the best fit values of each. Subsequently, 68% and 95% confidence level 
contours for various parameters have been presented, where all other parameters have been 
kept at their global best fit values (in contrast to most studies where they are fixed at the 
corresponding SM values). 

The most important conclusion we draw is that it is too early to say that signals reveal 
the standard model Higgs, since the 2a contours allow departure from the SM values. A few 
general trends that show up are 

• A fermiophobic Higgs is by and large disfavoured. 

• There is in general the hint of a relative sign between the couplings to the fermion and 
the gauge boson pairs. 

• A non-trivial phase in the top-quark coupling can have a rather important role. Inter- 
estingly, the SM values of the remaining parameters tend to get included within the 
2(7 contours once this phase is turned on. 

• It still seems possible to accommodate an invisible decay width of the Higgs. In par- 
ticular, it can be as large as 60% if one allows a breakdown of the custodial SU(2) 
symmetry. Though this is prima facie restricted by electroweak precision data, it re- 
mains to be checked whether the values of (3w 7^ l3z can be made consistent with 
precision data, by introducing more than one gauge invariant higher-dimensional op- 
erators at the same time. However, invisible branching ratios upto 20% or more are 
possible even without such manoeuvres. 

Data from the LHC are currently in a state of fiux, and therefore the numerical results 
on which our analysis is based can change with growing statistics. However, while the trends 
pointed out by us might change as more accurate data become available, the rather general 
approach used by us should continue to serve as a template for future analyses. 
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